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Abstract 

Thermal leptogenesis explains the observed matter-antimatter asymmetry of the uni- 
verse in terms of neutrino masses, consistent with neutrino oscillation experiments. 
We present a full quantum mechanical calculation of the generated lepton asymmetry 
based on Kadanoff-Baym equations. Origin of the asymmetry is the departure from 
equilibrium of the statistical propagator of the heavy Majorana neutrino, together 
with CP violating couplings. The lepton asymmetry is calculated directly in terms of 
Green's functions without referring to "number densities". Compared to Boltzmann 
and quantum Boltzmann equations, the crucial difference are memory effects, rapid 
oscillations much faster than the heavy neutrino equilibration time. These oscillations 
strongly suppress the generated lepton asymmetry, unless the standard model gauge 
interactions, which cause thermal damping, are properly taken into account. We find 
that these damping effects essentially compensate the enhancement due to quantum 
statistical factors, so that finally the conventional Boltzmann equations again provide 
rather accurate predictions for the lepton asymmetry. 



1 Introduction 



Standard thermal leptogenesis [T] provides a simple and elegant explanation of the origin 
of matter in the universe. Baryogenesis via leptogenesis naturally emerges in grand unified 
extensions of the Standard Model, which incorporate right-handed neutrinos and the see- 
saw mechanism, and the predicted connection between the cosmological matter-antimatter 
asymmetry and neutrino properties is in remarkable agreement with the present evidence 
for neutrino masses [2]. 

Leptogenesis is an out-of-equilibrium process in the high-temperature symmetric phase 
of the Standard Model. It makes use of nonperturbative properties of the Standard Model, 
the sphaleron processes which change baryon and lepton number [3], and it requires CP 
violation in the lepton sector and quantum interference in the thermal bath. Almost 
all quantitative studies of leptogenesis to date are based on Boltzmann's classical kinetic 
equations for the description of the nonequilibrium process [2j. 

In this article, we discuss a full quantum mechanical calculation of the generated lep- 
ton asymmetry based on Kadanoff-Baym equations [4J and the Schwinger-Keldysh formal- 
ism [HHT]. The main result has previously been reported in [8]. Here we give a detailed 
derivation of the result, discuss its interpretation and set the stage for future computations. 
Further work is still needed to obtain a 'quantum theory of leptogenesis' that can predict 
the cosmological matter-antimatter asymmetry in terms of neutrino properties without 
uncontrolable assumptions. 

Conventional leptogenesis calculations based on kinetic equations suffer from a basic 
conceptual problem: the Boltzmann equations are classical equations for the time evolution 
of phase space distribution functions; the involved collision terms, however, are obtained 
from zero-temperature S-matrix elements which involve quantum interferences. This is 
in contrast to other successful applications of the Boltzmann equations in cosmology, like 
primordial nucleosynthesis, decoupling of photons or freeze-out of weakly interacting dark 
matter particles, where the collision terms arise from tree-level S-matrix elements. In the 
case of leptogenesis, clearly a full quantum mechanical treatment is necessary to understand 
the range of validity of the Boltzmann equations and to determine the size of possible 
corrections [9j. 

In recent years, various attempts have been made to go beyond Boltzmann equations. 
In 12], a solution of Kadanoff-Baym equations for leptogenesis has been found to leading 
order in a derivative expansion in terms of distribution functions satisfying the Boltzmann 
equations. Various thermal corrections, in particular quantum statistical factors and ther- 
mal masses, have been included [101 - 113] . Quantum Boltzmann equations have been derived 
from Kadanoff-Baym equations for scalar and Yukawa theories [TUfTS] and for leptogene- 
sis [TBVfl^ . Except for |16) . they do not contain memory effects, but they yield the correct 
statistical factors which go beyond the Boltzmann equations |S1[T51[T71[IS1[2D]- Quantum 
Boltzmann equations have important applications for resonant leptogenesis [T6], flavoured 
leptogenesis |2T | |22 ] and A'"2-leptogenesis [23j . Similar techniques have been developed for 
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Figure 1: Evolution of heavy neutrino abundance Nn^ and lepton asymmetry Nb-l for 
typical leptogenesis parameter: Mi = 10^° GeV, fhi = 87iTi{vew/Miy = 10~^ eV, e = 
10~^; the inverse temperature z = Mi/T is the time variable. The dashed (full) lines 
correspond to thermal (vacuum) initial conditions for the heavy neutrino abundance; the 
dotted line represents the equilibrium abundance. From [30j . 



electroweak baryogenesis [23H27] and for coherent baryogenesis [28]. 

The quantum treatment of leptogenesis discussed in this paper is entirely based on 
Green's functions, thus avoiding all approximations needed to arrive at Boltzmann equa- 
tions. Our work is based on [29j, where the approach to thermal equilibrium has been 
discussed in terms of Green's functions for a toy model, a scalar field coupled to a large 
thermal bath. In leptogenesis it is the heavy neutrino which is weakly coupled to the stan- 
dard model plasma containing many degrees of freedom. The nonequilibrium propagator 
of the heavy neutrino is obtained by solving the Kadanoff-Baym equations. The induced 
quantum corrections of the lepton (and Higgs) propagators then yield the wanted lepton 
asymmetry. 

In general baryogenesis requires departure from thermal equilibrium. For the cosmo- 
logical baryon asymmetry, this is provided by the Hubble expansion of the universe and, 
possibly, also by initial conditions. This can be seen in Fig. [1] where the time evolution 
of heavy neutrino abundance and lepton asymmetry, as predicted by the Boltzmann equa- 
tions, are shown for two different initial conditions: thermal and zero heavy neutrino abun- 
dance. In the first case, the Hubble expansion leads to an excess of the neutrino abundance 
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at T ~ 0.3 Mi; shortly afterwards, washout processes are no longer in equilibrium and the 
lepton asymmetry is 'frozen in'. This is the standard out-of-equilibrium decay scenario of 
baryogenesis. In the second case, interactions with the thermal bath first bring the heavy 
neutrino into thermal equilibrium; due to the departure from thermal equilibrium during 
this time, an initial lepton asymmetry is generated. Around T ~ 0.3 Mi, this asymmetry is 
washed out and, as in the first case, the final lepton asymmetry is generated. Remarkably, 
the initial and the final asymmetry have about the same size. For the generation of the 
initial asymmetry the change of temperature due to the Hubble expansion is not impor- 
tant. This allows us to make a significant technical simplification in our analysis. Since 
our goal is the comparison of Boltzmann and Kadanoff-Baym equations, we concentrate on 
the computation of the initial asymmetry at constant temperature. We expect differences 
between the classical and the quantum approach to be of similar size in the generation of 
the final asymmetry. In our numerical analysis we shall consider temperatures T < M, 
where the heavy neutrino production rate is not strongly affected by the effect of thermal 
masses of lepton and Higgs fields |11H13] . 

We consider an extension of the Standard Model with additional gauge singlet fermions, 
i.e., right-handed neutrinos, whose masses and couplings are described by the Lagrangian 
(sum over 

^ = ^SM + l^iig'^Ri + lLi<PKj^Rj + 'l^jKjki4> - ^Mij {vR.jl'Rj + l^j^'ki) ■ (1-1) 

Here i/^ = Cz/^, C is the charge conjugation matrix and (p = ia24>*', SU(2) isospin indices 
have been omitted. For simplicity, we consider the case of hierarchical Majorana masses, 
^ Ml = M, and small Yukawa couplings of the lightest heavy neutrino A^i = A^, 
Ail ^ 1; such that the decay width is much smaller than the mass. Leptogenesis is then 
dominated by decays and inverse decays of A^, and it is convenient to integrate out the 
heavier neutrinos. From Eq. fll.ll) one then obtains the effective Lagrangian 

C =CsM + h^i^N + hi^X*,N + N^XnClLi(l) - ^MN^CN 

+ ^ViA<PClLj(p + ^v*ijki^CTL^<P , (1-2) 
with A^ = Ufii + u'^i, and the familiar dimension-5 coupling 

= 5^A,fc^Aj,. . (1.3) 

k>i ^ 

Using this effective Lagrangian has the advantage that vertex- and self-energy contributions 
to the CP asymmetry in the heavy neutrino decay |3T1433] are obtained from a single 
graph [9]. 

The paper is organized as follows. In Section 2 we present solutions of the Boltzmann 
equations for the heavy neutrino distribution function and the lepton asymmetry, which 
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are useful for later comparison with the KadanofT-Baym equations. Some results from 
nonequilibrium quantum field theory (QFT), in particular equilibrium correlation functions 
and Kadanoff-Baym equations, are recalled in Section 3. Section 4 contains some of the 
main results of this paper: analytic solutions of spectral function and statistical propagator 
for the heavy neutrino. These are needed for the computation of the lepton asymmetry, 
which is carried out in Section 5. A detailed comparison of the Boltzmann result and 
the KadanofT-Baym result is given in Section 6, and numerical results for the generated 
lepton asymmetries are compared in Section 7. Summary and conclusions are given in 
Section 8, and various details, including equilibrium correlation functions, Feynman rules, 
a discussion of the zero-width limit and the computation of some integrals are contained 
in Appendices A - D. 



2 Boltzmann equations 

The Boltzmann equations for the time evolution of the distribution functions of heavy 
neutrinos, lepton and Higgs doublets are well known [3l]. As discussed in the previous 
section, we focus on the generation of the 'initial asymmetry' (cf. Fig. [T]), which allows 
us to neglect Hubble expansion and washout terms and to work at constant temperature 
T. The distribution function of the heavy neutrinos is then determined by the first-order 
differential equatioij^ 



X [/;v(t,a;p)(l-/K^))(l + ^(?))-/K^)/<^(g)(l-/7v(t,^p))] , (2.1) 
with vacuum initial condition, 

/^(0,a;p)=0; (2.2) 



here Up = a/ JVP + p^, k and q are the energies of N, I and with equilibrium distribu- 
tion functions /; and f^, respectively; the averaged decay matrix element is \M{N{p) — )■ 
l{k)(j){q))\'^ = 2 (^X^Xj^^p ■ k (cf. [9]). For the momentum integrations we use the notation 

. (2.3) 



(27r)32a; ' ' ' ' 

In most leptogenesis calculations one directly computes the number density, 

nNit) = J -0^fNit,Up), (2.4) 



^To simplify notation, we use the same symbol for the modulus of 3-momentum and 4-momentum, e.g. 
k = |k| and k = (|k|,k). 
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assuming kinetic equilibrium. 

The sum of decay and inverse decay widths, whose inverse is the time needed to reach 
thermal equilibrium [33], is given by 

rp=(AU)^^— / {27r)'6\k + q-p)p-k fi^{k,q) , (2.5) 
where we have introduced the statistical factor (cf. [35j) 

q) = mfM + (1 - flikm + fM) 

= l-fi[k) + U{q) . (2.6) 

Neglecting the momentum dependence of the heavy neutrino width (Fp = F), one easily 
obtains the solution of the Boltzmann equation (12. ip with vacuum initial condition, 

/^(t,u;p) = /^«(a;p)(l-e-^*) , (2.7) 

where the equilibrium distribution is 



and /3 = 1/T is the inverse temperature. 

To compute the lepton asymmetry, we need the Boltzmann equation for the lepton 
distribution function, 



{2'RY6\k + q-p) 
X [|M(/0 ^ N)\^mU{q){'^ - fNit,co^)) 

- \M{N^l<P)\'Mt,u^){l-fi{k)){l + U{q))] , (2.9) 

where now (9(A^) corrections to the matrix elements have to be kept. Using Eq. (12. 7p one 
obtains for the lepton asymmetry 

fLi{t,k) = Ut,k)-fu{t,k) , (2.10) 

with initial condition /Li(0, k) = 0, 

fu{t, k) = -eul [ {2T:Y5\k + q-p)p.k f,^{k, q)r^'{oo^)^ (l - e"^*) , (2.11) 



where we have defined 

3 
16n 



e., = — lm{A*i(r^A*),i}M . (2.12) 
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Summing over all lepton flavours, the generated lepton asymmetry is proportional to the 
familiar CP asymmetry [9], 



4- (AtA),, 16vr (AtA),, " ^ • ^ 

For later comparison with solutions of the Kadanoff-Baym equations, it is convenient 
to rewrite Eq. (12.1 1|) as a 4- fold integral, 

Mt, k) = -eu^-^ [ k-k' {2nY6\k + q - p){2nY5\k' + q' - p) 

x/.,(A:,g)/^>p)^(l-e-^*) . (2.14) 

The integrand is now proportional to the averaged matrix element \M{l(f) — )■ /0)P = 
2k ■ k'{X^X)ii/M^ (cf. [9j), which involves the product of the 4-vectors k and k'. At low 
temperatures, T ^ M, the integrand falls off like e~'^'^p < e~^*^, i.e., the generated 
asymmetry is strongly suppressed. In standard leptogenesis calculations one considers the 
integrated lepton asymmetry, 

--=Ey (^/-(^'^)- (2-15) 

The number densities (12. 4p and Ul correspond to the comoving number densities Nn^ 
and |Nb-l| shown in Fig. [H in the initial phase of the time evolution, i.e., for T > 0.3 M. 



3 Nonequilibrium QFT and Kadanoff-Baym equations 

In the following, we briefly introduce concepts and quantities from nonequilibrium quantum 
field theory that are necessary for our computation (cf. [5B|IB7] ). A thermodynamical system 
is represented by a statistical ensemble described by a density matrix g. The expectation 
value for an operator A is then given by 



{A) = TT{gA) , (3.1) 

where we have adopted the usual normalisation Ttq = 1. Solving the initial value problem 
for g allows to compute all observables for all times. Direct computation of the time 
evolution of g is difficult. Generically, the von Neumann (or quantum Liouville) equation 
of motion for g can only be solved perturbatively for a reduced density matrix with an 
effective Hamiltonian. In most practical applications to date, a number of additional 
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assumptions are made that lead to effective Boltzmann equations, which can take account 
of coherent oscillation^, or quantum corrected Boltzmann equations (cf. Section 6)|f|. 

Instead of the time evolution of the density matrix, one can also directly study the 
equations of motion of the correlation functions of the theory. The infinitely many degrees 
of freedom of the initial density matrix are then mapped onto their infinitely many initial 
conditions. Though a full characterisation of the system in principle involves all n-point 
functions, it is often sufficient to study the one- and two-point function. This applies to 
the problem considered in this work. 



3.1 Correlation functions for lepton and Higgs fields 

Leptogenesis occurs at temperatures above the electroweak scale where sphaleron processes 
are active and transfer the generated lepton asymmetry to a baryon asymmetry. Hence, 
the Standard Model is in the symmetric phase and the four real degrees of freedom of the 
Higgs doublet correspond to four massless real scalar fields. 

The spectral function and statistical propagator of a real scalar field 0, A~ and A"*", 
respectively, are defined as 

^-{xi,X2)=imxi)A{x2)]) , (3.2) 

A+(X1,X2) = ^({0(X1),0(X2)}) . (3.3) 

Here only contributions from connected diagrams are to be included to compute the dressed 
correlation functions. These fulfill the symmetry relations 

A-(xi,a;2) = -A-(a;2,xi) , (3.4) 
A+(xi,a;2) = A+(x2,a;i) , (3.5) 

which follow directly from the definitions. 

The functions A^ have an intuitive physical interpretation. The spectral function A~ 
is the Fourier transform of the spectral density, 

P^{t,uJ) = -^ j ge^-^A-(t+|,t-|), (3.6) 

where we have used the relative and total time coordinates, y = ti — t2 and t = (ti + 12)/2, 
respectively. 

The spectral density pq{t,uj) characterises the density of quantum mechanical states in 
phase space. Propagating states, or resonances, appear as peaks in the spectral function. 



^See [38. 39] for an application to neutrino oscillations. 

■^In [40P4T] an approach based on first principles has been suggested that is applicable if the occupation 



numbers for the out-of-cquilibrium fields are small. 
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The statistical propagator contains the information about the occupation number of each 
state. 

In the following we shall also need the Wightman functions 

A>(Xi,X2) = (0(X1)0(X2)) , (3.7) 
A<(Xi,X2) = (0(x2)0(xi)) , (3.8) 

which are related to by 

A"(xi,X2) = i (A>(xi,X2) - A<(xi,X2)) , (3.9) 

A+(xi,X2) = ^(A>(xi,X2) + A<(xi,X2)) . (3.10) 

Using microcausality and the condition for canonical quantization, 

[0(Xi),0(x2)]|t,=t, = [0(Xi),0(x2)]|t,=t, = , (3.11) 
[(l){Xi),<p{x2)]\t^=t2 =i5(xi-X2) , (3.12) 

one obtains boundary conditions in y = ti — 12 for A~, 

A-{xi,X2)\t,=t,=0 , (3.13) 
9t,A^(xi,X2)|ti=t2 = -dt2A'{xi,X2)\ti=t2 = 5(xi -X2) , (3.14) 
dt,dt2A-{xi,X2)\ti=t2 = . (3.15) 

Note that these conditions do not depend on the physical initial conditions of the system 
encoded in the initial density matrix. These enter via the initial conditions for the statistical 
propagator. 

Analogous to A='=, one can define the spectral functions and statistical propagators for 
fermions. The fermionic fields in the Lagrangian (11. 2p are massless left-handed leptons 
(Weyl fields lu) and a massive neutrino (Majorana field A^). For the massless leptons, 
spectral function and statistical propagator are defined as 

iSLij)a(3{Xl,X2) = i{{lLia{Xl)jLjl3ix2)}) , (3.16) 
iStij)ap{Xl,X2) = ^{[lLia{Xl)jLjp{x2)]) , (3.17) 

where a and /3 are spinor indices, and SU(2) indices were omitted for notational simplicity. 
The subscript L denotes the projection to left-handed fields, i.e., = PlS*^, where 
Pl = {1 — 7^)/2 and 5*^ are the propagators for Dirac fermions. As for bosons, we shall 
need the functions 

iSLij)a(3iXi,X2) = {lLia{Xl)lLjp{x2)) , (3.18) 
{Stij)al3{Xl,X2) = -(jLjl}{x2)lLia{Xl)) , (3.19) 



which are related to spectral function and statistical propagator by 

^Liji^i, = i {Sujixi, X2) - S^ij{xi, X2)) , (3.20) 
Stijixu X2) = I {Sl^ixu X2) + S<j{x,, X2)) . (3.21) 

The propagators have the symmetry properties 

7o [SZijixi, X2)] ^ 7o = -SZji{x2, xi) , (3.22) 
70 [Slj{x^,X2)] ^ 70 = St.,{x2, X,) . (3.23) 

The canonical quantization condition, 

{kiaixi), l[jgix2)} = PlafiSijSixi - Xg) , (3.24) 

implies the boundary condition for the spectral function 

S'Eij(xi,X2)\ti=t2 = iPLSij5(-xi - X2) . (3.25) 
Finally, spectral function and statistical propagator for the Majorana field N read 

G-p{xuX2) = t{{N,,{x,),N0{x2)}) , (3.26) 
GU^i,X2) = ^{mx,),N^{x2)]) . (3.27) 

They have the symmetries 

G-{xi,X2) = G-{x2,xi f , (3.28) 
G+{xi, X2) = -G+{x2,x,f . (3.29) 

The canonical quantization condition, together with the Majorana property N = CN'^, 
implies the boundary condition 

G-{xi,X2)\t,=t, = ^7°5(xl - X2)C-i . (3.30) 

As for scalars, the physical initial conditions enter as boundary conditions for the statistical 
propagator. In the following, we will consider two types of initial conditions, thermal 
equilibrium and Gaussian initial correlations, for which we solve the equations of motion 
in the following section. Analogous to real scalars, the functions are defined as 



G>p{xi,X2) = {N^{xi)N^{x2)) , (3.31) 

G<^{xi,X2) = -{Np{x2)N^{xi)), (3.32) 

with the usual relations to spectral function and statistical propagator, 

G-ix,,X2)^iiG>{xi,X2)-G<{xi,X2)) , (3.33) 

G+{x,,X2) = ^{G>{xuX2) + G<{x,,X2)) . (3.34) 
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Figure 2: Path in the complex time plane for nonequilibrium Green's functions. The 
contour runs from some initial time x° = ti + ie {U = 0) parallel to the real axis {x^ = t + ie) 
up to some final time tf+ie and returns to ti—ie. To compute physical correlation functions 
for arbitrary times t > ti, one takes the limits — )■ oo and e — > 0. 

3.2 Equations of motion 

In thermal leptogenesis, the deviation from thermal equilibrium that is necessary to create 
a matter-antimatter asymmetry is due to the heavy Majorana neutrinos which are out of 
equilibrium. The equations of motion for their correlation functions can be obtained 
via the Schwinger-Keldysh formalism [6]. The basic quantity is the Green's function with 
time arguments defined on a contour C in the complex x°-plane, known as the Keldysh 
contour (cf. FigJ2]), 

Gcixi, X2) = edxl x°)G'>(xi, X2) + edxl x?)G'<(xi, x^) . (3.35) 

Here the ^-functions enforce path ordering along the contour C. The necessity of considering 
Green's functions with time arguments on the Keldysh contour (rather than the real axis) 
is a consequence of the fact that nonequilibrium processes are initial value problems. The 
system is prepared at initial time ti, its state at later times is unknown. Hence, the usual 
approach to define a S-matrix by projection onto asymptotic 'in' and 'out' states, sending 
initial and final time to infinity, cannot be applied. When using the Keldysh contour which 
starts and ends at the same time tj^, no knowledge of the system's state at t = ±00 is 
needed to define a generating functional for correlation functions. 
The Green's function Gq satisfies the Schwinger-Dyson equation 

C {igli - M)G e{xi,X2) - i j d'^x'CT,cixi,x')Gcix',X2) = iSdxi - X2) , (3.36) 

where CSc(xi,x') is the self-energj|f| on the contour and ^1 = 'j'^d/dxi- Like the Green's 
function, also the self-energy can be decomposed as 

Jlc{xi,X2) = Oc{xl,X2)Tj>{Xi,X2) + 6c{x2,X°^)Tj^{Xi,X2) . (3.37) 

'^Due to this fact this formahsm is sometimes called 'in-in' formalism, in contrast to the 'in-out' for- 
malism used to compute the S-matrix. 

^An explicit factor C is factorized for later convenience. 
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In the Schwinger-Dyson equation (13.36^ the time coordinates of Gq and Sc can he on the 
upper or the lower branch of the contour. 

The familar time-ordered Feynman propagator is obtained from Gc{xi,X2) when both 
time arguments he on the upper branch, and therefore denoted by G^^. Correspondingly, 
Gc{xi,X2) with both time arguments on the lower part of the contour corresponds to an 
anti-time-ordered propagator, denoted as G^^. For correlators with one time argument on 
the upper and one on the lower part of the contour, referred to as G^^ and G^^, the order 
of field operators is fixed by the path ordering: operators on the upper branch are always 
'earlier' than those on the lower branch (cf. I3.35p . Altogether, one has 

G^2(xi,X2) = G<(xi,X2) , (3.38) 
G21(xi,X2) = G>(xi,X2) , (3.39) 

G^^(xi,X2) = G+(xi,X2) - -sign(x° - X2)G"(xi,X2) , (3.40) 

G22(xi,a;2) = G+(xi,X2) + ^sign(x?-a;°)G-(xi,X2) ; (3.41) 

the last two relations are easily verified by inserting the definitions of G^. 

In a perturbative expansion of the Schwinger-Dyson equation ( I3.36P in terms of Feyn- 
man diagrams, time arguments of internal vertices can lie on either branch. Hence, the 
number of contributing graphs doubles with each internal vertex since this can lie on the 
upper or the lower branclu. Two upper vertices are connected by G^^, two lower vertices by 
G^^ and vertices of different type by G^^ and G^^. Each lower vertex leads to an additional 
factor —1. 

Like the Green's function, also the self-energy Sc, the sum of all one-particle irreducible 
graphs, can be dissected into components Tl'\ with k and / being 'contour indices' as defined 
above. Analogous to fl3.38p and fl3.39p one then defines self-energies and, following 
(I3.33P and fl3.34p . self-energies via the equations 

S"(xi,X2) = i (S^(xi,X2) - S'^(xi,a;2)) , (3.42) 

S+(xi,X2) = ^(S>(xi,X2) + S<(xi,X2)) . (3.43) 

Since the self-energies S'^' are directly related to the full Green's functions G^\ they also 
satisfy the relations (lOHl) - ([331]). 

Using the above relations for G^^ and E'^', one obtains, after a straightforward cal- 
culation, from the Schwinger-Dyson equation (I3.36P a system of two coupled differential 
equations for Gp, the Kadanoff-Baym equations. Due to spatial homogeneity, we can 
consider the equations for each Fourier mode separately, 

G(«7°5t, -p7-M)Gp(ti,t2) =- r dt'CT.^{h,t')G^{t\h) , (3.44) 



^This fact is sometimes referred to as 'doubling of degrees of freedom'. 
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C(z7°5t, - P7 - M)G+(ti,t2) = - / dt'CJ:;ih,t')G^it',h) 

Jti 

+ / c/t'CSp(ti,t')G'+(t',t2) • (3.45) 

For the lepton propagators S^y_ one obtains the same equations, with CSp replaced by the 
lepton self-energies 11^ and no charge conjugation matrix C multiplying the kinetic term. 

The Kadanoff-Baym equations f l3.44p and (I3.45P are exact. They contain all quantum 
and non-Markovian effects including the dependence on the initial time tj. Furthermore, 
in contrast to usual linear response techniques, they do not rely on any assumption re- 
garding the size of the initial deviation from equilibrium. The equations in this form are 
valid for arbitrary nonequilibrium initial states which can be parameterized by Gaussian 
initial correlations. This covers the case considered in this work since the generated lepton 
asymmetry involves to leading order in the Yukawa coupling only the 2-point functions 
of the heavy neutrino. When higher order initial correlations play a significant role, the 
Kadanoff-Baym formalism is still applicable, but the equation for the statistical propagator 
contains extra terms at tj [42j. In [36j, thermalization has been studied for a scalar field 
theory using the equation of motion for the statistical propagator. 

In nonequilibrium quantum field theory, instead of distribution functions, quantum 
mechanical correlation functions characterise the state of the system. The interactions 
enter via the self-energies which, via the generalized cutting rules, contain all possible 
processes. Encoding this information in the self-energies avoids potential problems related 
to the definition of asymptotic states for unstable particles as well as the substraction 
of real intermediate state contributions in Boltzmann equations. Note, finally, that the 
integro-differential equations (I3.44p . f l3.45p do not suffer from the late time uncertainties 
or secular terms that perturbative expansions of Boltzmann equations are often plagued 
with when applied to multiscale problems (cf. [36j). 

3.3 Weak coupling to a thermal bath 

The Kadanoff-Baym equations provide a tool to study the dynamics of arbitrary nonequi- 
librium systems. Unfortunately, in most cases they can only be solved numerically. As 
discussed in the introduction, in this work we consider a rather simple system: one field 
that is out of equilibrium (A^) is weakly coupled to a large thermal bath of Standard Model 
fields. This leads to a number of simplifications compared to the general case that allow to 
find analytic solutions. We have previously studied scalar field models of this type [211113]. 
Here we extend the methods developed therein to the case of thermal leptogenesis. 

The Standard Model interactions keep the bath in thermal equilibrium. The corre- 
sponding time scale Tsm ~ ^/{q'^T) at temperature T ~ M is much shorter than the 
equilibration time r^r ~ 1/ (A^M) of the heavy neutrino, which governs the generation of 
the lepton asymmetry: Tsm ^ 'Tn- Lepton number changing processes in the thermal bath 
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Figure 3: One- and two-loop contributions to the lepton self-energy correspoding to 
washout terms, a) - d), and CP violating terms which generate a lepton asymmetry, e) 
and f). 

are shown in Fig. |3l As in the case of Boltzmann equations discussed in Section 2, we 
focus on the CP violating interaction generating the lepton asymmetry that correspond to 
Fig.Se) and f). 

To evaluate these graphs we need the correlation functions of lepton and Higgs fields 
in the thermal bath. A system in thermal equilibrium is described by the density matrix 



where "H is the Hamiltonian of the system, /3 is the inverse temperature, Qi are conserved 
charges and /ij are the corresponding chemical potentials. As expected for an initial state 
after inflation, we set all chemical potentials equal to zero. 

Equilibrium correlation functions of a spatially homogeneous system only depend on 
space-time differences, and it is convenient to consider the Fourier transforms. 



Qeq 



exp if3 {-n + fiiQi)) 



(3.46) 




(3.47) 



(3.48) 
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Figure 4: Path in the complex time plane for equilibrium correlation functions. 



The equilibrium density matrix fl3.46p then corresponds to a shift in imaginary time. This 
leads to the well-known Kubo-Martin-Schwinger (KMS) relations (cf. |44) ) 

A<(a;) = e-^-Aliuj) , S<{uj) = -e-^-S>{u;) , (3.49) 

which imply 

A+M = + /^m) AqM = -^coth A^M , (3.50) 

S^ioo) = - M^)) S^ioo) = tanh (^^) 5, (c) , (3.51) 

where 

/^M = ;^' /'H = ;^^ (3-52) 

are Bose-Einstein and Fermi-Dirac distribution functions, respectively. Note that the en- 
ergy uj is not on-shell. 

Equilibrium Green's functions can be calculated in the real-time formalism using the 
contour in the complex time plane, which is shown in Fig. |H For the free equilibrium 
propagators of massless lepton and Higgs fields one obtains (g = |q|, k = |k|, cf. [44j ) . 

A^iy) = -sm{qy) , (3.53) 
q 

K(y) = Yq "^^^^ ( y) "^^'^^^^ ' ^^-^^^ 

Sl\,{y) = Pl {ilo cos(%) - ^ sm{ky)^ , (3.55) 
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Figure 5: One-loop contribution to the self-energies CSp of the Majorana neutrino A^. 

All other propagators can be obtained as linear combinations using the relations described 
in the previous paragraph. A complete list is given in Appendix A. 

In the following sections we shall see that the calculation of the lepton asymmetry 
represents an initial value problem which can be treated based on the real time formalism 
together with the Keldysh contour Fig. O Thermal and nonthermal properties of the 
system are then encoded in the initial values of the various Green's functions. 

4 Nonequilibrium correlation functions 

The assumption of weak coupling to a large thermal bath with negligible backreaction in the 
framework of Kadanoff-Baym equations implies that self-energies for the heavy neutrinos 

are computed from equilibrium propagators of bath fields only. This also corresponds 
to a leading order perturbative expansion in the coupling constant. 

Perturbative expansions of Boltzmann equations in multiscale problems are known 
to suffer from uncertainties, so-called secular terms, at late times. The Kadanoff-Baym 
equations fl3.44p and fl3.45p in full generality are free of secular terms and consistently 
include all memory effects. Nevertheless, the neglect of backreaction in the computation 
of S corresponds to a truncation in the perturbative expansion in the Yukawa couplings A, 
which might introduce similar uncertainties related to the multiscale nature of the problem. 
However, in the system of consideration contributions of higher order in A are not only 
suppressed by the smallness of the coupling, but also by the number of degrees of freedom 
in the bath that justify the neglect of backreaction. Hence, we expect potential problems 
due to secular terms not to be relevant. 

The assumption that the background medium equilibrates instantaneously on the time 
scale of the asymmetry generation leaves open the details of the equilibration process. In 
reality, there are effects related to the finite equilibration time and the finite size of the 
quasi-particles. As we shall see in Section 5, these quantities play a crucial role in the 
Kadanoff-Baym result for the lepton asymmetry. 

The self-energy for the heavy neutrino N to leading order in A is given by the diagram 
in Fig. [5l It contains time-translation invariant propagators of bath fields only, and hence 
it is also time-translation invariant. As shown in [29j, this implies that also the spectral 
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function is time-translation invariant, Gp{ti,t2) = Gp{y), y = ti — t2. In this case we can 
find the general solutions to the Kadanoff-Baym equations without further approximations. 



4.1 Equation for the spectral function 

Let us now consider the equation for the spectral function of the Majorana neutrino. After 
an obvious change of variables, the Kadanoff-Baym equation (I3.44p becomes, 

C{i^'dy-^-M)G-{y)- r dy'CE-{y-y')G-{y')=0 . (4.1) 

Jo 

Defining the Laplace transform 



Gp(s)=/ dye-^yG^{y), ^^(s) = dye-^yT.^{y) , (4.2) 
Jo Jo 

one obtains from Eq. (14. ip 

(^,/, - irr - M - Sp (s)) Gp (s) = ^T^Gp (0) . (4.3) 

Using the boundary condition (I3.30p . 

G-{Q)=i^'G-\ (4.4) 



p 

this leads to 

G'-(s) = - (^i7°s-p7-M-S-(s)) ^ G-^ . (4.5) 
The inverse Laplace transform is given by 

Gp(?/) = ^ |^e^^Gp(.), (4.6) 

where Cb is the Bromwich contour (see Fig. [6]): The part parallel to the imaginary axis 
is chosen such that all singularities of the integrand are to its left; the second part is the 
semicircle at infinity which closes the contour at Re(s) < 0. 

From the definition of the Laplace transform one can see that the self-energy Sp(s) is 
analytic on the real s axis, but has a discontinuity across the imaginary axis. This gives 
rise to the spectral representation 

. dp, S-(po) 

Sp(^5) = ^ / 1^— • (4-7) 

Note that the retarded and advanced self-energies are given by 

Sp(-^u; + 6) = E» , (4.8) 
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Figure 6: Bromwich contour 



These self-energies are determined by the discontinuity of Sp [s] 



discSp {—iuj) 



Sp(-iu; - e) 



with the real part given by the principal value, i.e., 



S„ (—iu ± e) = iV 



(4.9) 



(4.10) 



(4.1i; 



This representation of the self-energy is familiar from the theory at zero temperature. 

We are now ready to calculate the spectral function in terms of the self-energy Sp(ci;). 
Its Laplace transform has singularities only on the imaginary axis. Hence the Bromwich 
contour can be deformed as Cb — + fiJ^^^ (see Fig. [6]), which yields for the spectral 
function 



^^e(-+^)^Gp(za; 



-oo 
f oo 



■e)+ r°°^e(— )^Gp(.c.-e) 



— e-^(Gp(-.a; + e)-Gp(-^a;-e)) . 



(4.12) 



The Fourier transform of the spectral function, 

POO 

p^iu) = -I / dye'^yG-^{y) , 



(4.13) 
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is then given by 



\ 



—I 



(4.14) 



Here we have assumed that the divergent contribution of the real part has aheady been 
absorbed into mass and wave function renormahzation, so that Pp{uj) represents the renor- 
mahzed spectral density. The finite part of the self-energy is negligable because of the 
small Yukawa coupling. 

A straightforward calculation yields for the self-energy (cf. [33]). 



(4.15) 



where we have defined 

ct{p; k, q) = fi4k, g)(27r)^ {6\p - k - q) + 5\p + k + q)) 

+ Mk, q)(2TxY {5\p + k~q) + 5\p -k + q)) 

with the statistical factors 



fi^{k, q) = l- fi{k) + f^{q) , f^{k, q) 



+ m . 



(4.16) 



(4.17) 



Note that k and q are on-shell, i.e., k = {k, k) and q = {q, q), whereas p = {u, p) is off-shell. 
The properties of the Dirac matrices and rotational invariance imply 



Sp (uj) = iap{u)-f° + z6p(a;)p7 , 



where 



ip(u;) =2(AU)^^ / ka{p;k,q) 

Jq,k 



6p(a;) = -2 (AtA)^^ 



1 



P Jq,k 

These functions satisfy the relations 



pk a{p; k, q) 



ap(wp) , hp{-Up) = -bp{up) . 



(4.18) 

(4.19) 
(4.20) 

(4.21) 



Using Eq. (I4.18P and linearising the denominators in Eq. fl4.14p in the small quantities 
ap{u) and bp{u), one obtains for the spectral density 



'+M) C-^ 



(4.22) 
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where 

uTp{u) = uap{u) + p^6p(w) 



= 2(AU)^ / p-ka{p;k,q) . (4.23) 

Jq,k 

On-shell, only the first of the 5-functions in a{p] k, q) contributes, and one obtains the 
width appearing in the Boltzmann equations, 

rpK)= (AtA)^^— / p.kfi^{k,q){27r)'6\p-k-q)^Tp, (4.24) 

which satisfies the relations 

rp(-a;p) = r_p(a;p) = rp(a;p) . (4.25) 

In the zero-width limit the spectral function fl4.22p reduces to the familiar expression in 
vacuum, 

Pp{uj) = 27rsign(a;)5(p2 - M'^){^+ M)C-^ . (4.26) 

The spectral propagator is now obtained by evaluating the Fourier transform of the 
spectral function fl4.22p . 

Gp{y)=^ r^e-Vp(u;), (4.27) 

which yields the final result 

Gp{y) = [iioOO^M + e-^.\y\nC-' . (4.28) 

Compared to the free spectral function only an exponential damping factor appears. This 
is a feature of the narrow-width approximation, analogous to the scalar field case discussed 
in m. 



4.2 Equation for the statistical propagator 

We now proceed to the solution of the second Kadanoff-Baym equation fl3.45p which, 
choosing = 0, reads 

fti 

C(^7°5t, - p7 - M)G;ih,t2) - / rft'CSp(ti - t')G;it',h) = Cp(ti - t2) , (4.29) 

^0 
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with the source term 

-12 



Cpih - t2) = - / ' dt'J:;{ti - t')G^{t' - t2) . (4.30) 

The general solution of f l4.29p takes the form 

Gp{ti,t2) = Gp{ti,t2) + Gp^^^^(ti,t2) , (4-31) 
where Gp{ti,t2) is the general solution of the homogeneous equation 

C{t-f%, - p7 - M)G;{h,t2) - r dt'GE-ih - t')G;{t',t2) = , (4.32) 

Jo 

and the 'memory integral', which contains non-Markovian effects, is given by 

G;„{ti,t2) = r dt' r dt"Gp{h - miit' - t")Gp{t" - t2) . (4.33) 
Jo Jo 

One easily verifies that the memory integral is a special solution of the inhomogeneous 
equation. 

In order to evaluate the memory integral we perform a Fourier transform of the self- 
energy {y = ti-t2), 

^p,mem(^l5^2) = 

/ ^ {£ dy^Gpiyiy^'-) Ki^) dy2Gp{~y2)e-^-'^^ e^^ . (4.34) 

Since the self-energy is computed with fields in thermal equilibrium, it satisfies the KMS 
condition (cf. (I3.5ip ) 

S+(a;) = -ltanhf^')s-(a;). (4.35) 



2 V 2 ^ 

Using the expressions fl4.18p and f l4.28p for self-energy and spectral function, respec- 
tively, which were derived in the previous section, it is now straightforward to calculate 
the memory integral explicitly. Neglecting terms OiTp) in the numerator, one finds 



27° [(wp sin(u;pt) + cos(wpt))e*(^+*^p/2)* - too] 
+ [tio ( sin(a;pt) - Up cos(a;pt)) e'("+*^''/2)t ^ \ ^-i ^ 

Up I 
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dye-'G;(-y) = _ _ .^p^^,,. x (4.37) 

[ico ( sm(a;pt) - cUp cos(a;pt)) e'^'^+'^p/')* + Wp] ) C"^ . 



M - p7 

+ — 



Up 



After inserting these expressions in Eq. (14.341) one can perform the w-integration using 
Cauchy's theorem. The integrand has two poleCl in the upper-half plane at a; = iTp/2±ujp, 
and two poles in the lower-half plane uj = —irp/2±ujp. The choice of the contour depends 
on the sign of the time variables in the exponent. The result is a sum of the contributions 
from all four poles. The expressions appearing in the numerator can be simplified by means 
of Eqs. f l4.2ip and f l4.24p for self-energy and equilibration width, respectively, 



(,o ^ ^).,,-.,(/ _ ^) - .r,(/ - ^£^) . (4.3.) 



Using these expressions one finally obtains for the memory integral, changing variables 
from (ti,t2) to (t,y), 

Gt^.Jt,y)= (4.40) 
- itanh (^loMupy) - ^^-^ cos{upy)^ (e-^'^l^l/^ - e'^'^') C'^ . 

Asymptotically, for t — t- oo, the memory integral becomes 

Gp^{t,y) = -^tanh (^^) (.TosinKy) - ^^^^ cosKy)) e~''^\y\/'C-' . 

(4.41) 

One easily verifies that Gp'^'^{t,y) indeed represents the equilibrium statistical propagator. 
For the Fourier transform one obtains 



/oo 
dye'-yC^'^^y) 
-oo 



1 ^ ^ f I3u\ 2wrp(a;) 



- tanh — — ^ (]/+ M) C- 

Uanh ppico) (4.42) 



^There are further poles at ujn = ±i7r(l + 2n)//3, n integer. However, their contribution to Gp^^^^ is 
0{Tp/M) and therefore neghgible. 
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= -|tanh(^)GpM, (4.43) 

i.e., the KMS condition (cf. fl3.5ip ) is indeed satisfied. 

In order to obtain the general solution of the inhomogeneous Kadanoff-Baym equation 
we have to add to the memory integral the general solution of the homogeneous equation 
( I4.32p . This equation is identical to the Kadanoff-Baym equation for the spectral function 
(14. ip with t2 playing the role of an additional parameter. Hence, the functional dependence 
of Gp{ti,t2) on the first argument ti can be obtained in the same way as for the spectral 
function. Applying the Laplace transform to (I4.32p one finds 

G+(.,t2) = : ^— -^7oG+(0,t2) . (4.44) 

The inverse Laplace transform then gives 

^+(^1,^2) = -G'p(ti)«7oG+(0,t2) . (4.45) 

The function Gp(0,t2) can now be determined by the symmetries f l3.28p and f l3.29p of 
Gp(ti,t2) , which imply 

G;{tut2f = -G;it2,h) . (4.46) 
This yields the result 

G+(ti,t2) = -Gp(ti)C7oG+(0,0)7oC^-'Gp(-t2) , (4.47) 

where Gp{0,0) is an antisymmetric matrix. 

Let us first consider the case of thermal initial condition, 

0(0, 0) = tanh ( ^] C-\ (4.48) 



2ujp V 2 



From Eq. (I4.47P one then obtains 

GT(ti,t2) = -I (iloMu^^y) -^^^cosiupy)^ tanh (^^) e-r<^(*^+*^)/^C-i . 

(4.49) 

Adding this expression to the memory integral Gp j^^j^ one obtains the equilibrium statis- 
tical propagator G^'^^ which is independent of t = {ti + t2)/2. Hence, as expected, the 
equilibrium statistical propagator is a solution of the full Kadanoff-Baym equation. 
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We are particularly interested in the case of vacuum initial condition, which corre- 
sponds to zero initial abundance for heavy neutrinos in the Boltzmann case. The vacuum 
propagators are obtained from the equilibrium ones in the limit /3 — > oo. Hence we choose 

^+^-(0, 0) = ^~J^ C-^ ■ (4-50) 

From Eqs. (14.3ip . (I4.40p and f l4.47p one then obtains the full solution for the statistical 
propagator, which interpolates between vacuum at t = and equilibrium for t — > oo. 



'itanh(^) e-^pl^l/^ + Z-We- 



C-\ (4.51) 



This result will be the basis for the calculation of the lepton asymmetry in the next section. 
All heavy neutrino propagators can be obtained as linear combinations of the spectral 
function G~{y) and the statistical propagator Gp{t,y). A full list is given in Appendix A. 

Finally, let us emphasize that the solution of the Kadanoff-Baym equation for the 
statistical propagator is not related to the equilibrium propagator by a simple change of 
the distribution function from /^(w) to some nonequilibrium function f"f^{t,u). This is in 
contrast to the assumption made in the derivation of Quantum Boltzman equations |16 [ [T7 t 
[T9l|21]. For a system close to equilibrium this assumption leads to a valid approximation 
of the Kadanoff-Baym equations [9j, but in general it is not justified. 

5 Lepton asymmetries 

We are now ready to calculate the lepton asymmetry which is generated during the ap- 
proach of the heavy Majorana neutrino N to thermal equilibrium. Our starting point is 
the flavour non-diagonal lepton current, which is obtained from the statistical propagator, 

ji'^{x) = -tT[YSt,^{x,x%,^,. (5.1) 

Since we consider a spatially homogeneous system, S^j{x, x') only depends on the difference 
X — x', and it is convenient to perform a Fourier transform. The zeroth component of the 
current, the 'lepton number matrix', is given by 

Lk.,(t,t') = -tr[7o5+k.,(t,0] • (5.2) 
One easily verifies that for free fields in equilibrium 

L^u{t,t) = fu{k) - h{k) , (5.3) 
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^3 

Figure 7: Two- loop contributions to the lepton self-energies 11^, which lead to a nonzero 
lepton number densities. 



where fn and ff^ are the distribution functions of leptons and anti-leptons, respectively. 

The lepton number matrix L]f_ij{t, t') can be directly computed from the self-energy cor- 
rections to the statistical propagator shown in Fig. |3l the external lepton couples to Majo- 
rana neutrino and Higgs boson, and also to Higgs boson and Higgs-lepton pair. Complex 
Yukawa couplings and quantum interference then lead to a non- vanishing lepton asymme- 
try. 

For a homogeneous system, the Kadanoff-Baym equation for the statistical propagator 
(cf. fl3.45p ) yields for each Fourier mode the equations 

- k7)5+k(t,t') = f dU\]^{t,t,)Sl^{U,t') 
Jo 

t' 

dt,UUt,h)SZ^ih,t') , (5.4) 
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St^it,t')i-^^'^^, -k7) 



dhS+^{t,h)T]-{tut') 



+ [ dtiSL^,it,ti)T]^{ti,t') . 
Jo 



(5.5) 



One then obtains for the time derivative of the lepton number matrix, dropping flavour 
indices (cf. [18j)H 

dMt, t) = lii [(«7o9t + l-1MSt^.it. t')]t=t' 



iii 
iti 



{i^,dt - k7)S+k(t,t') + SUt.mi'd, +k7) 



t=t' 



[ dt^W^{tM)St^{h,t')- f dt,T]+{t,t,)S^^{t,,t') 

Jo Jo 

+ [ dtiS^^^it,t,)T]^^{t,,t') - [ dt^S^^{t,t,)U+{h,t') 

Jo Jo 



t=t' 



Using properties of the trace and the identity between integration domains 



ti 



dti dt2 h / dt 



t2 



dti 



dti / dt 



one finds 



L^{t,t)=t[ dti [ dt2 tr[nk(ti,t2)5+k(t2,ti)-n+(ti,t2)5^k(i2,ti 
Jo Jo 



(5.6) 



(5.7) 



(5.8) 



Note that TT^ and are self-energies and propagators of the full theory including gauge 
interactions of lepton and Higgs fields. 

Using the relations for propagators and self-energies 



•^Lk — 2 ^^^^ ~^ ^^^^ ' ~ ^ ^^^^ ^^^^ ' 



(5.9) 
(5.10) 



one obtains from Eq. (15. 8 p an equivalent useful expression for the lepton number matrix, 

L^{t,t) = -! dti f c/t2tr[n>(ti,t2)5^k(^2,ti) -n<(ti,t2)5^k(^2,ti)] • (5.11) 
Jo Jo 

We want to calculate the lepton asymmetry to leading order in the small Yukawa cou- 
pling A, which can be achieved in a perturbative expansion. For the heavy neutrino propa- 
gator appearing in the loop, the departure from the equilibrium propagator is important]^ 



^We thank C. Weniger for helpful discussions. 

^We show in Appendix D that the equilibrium part of the propagator does indeed not contribute to 
the asymmetry. 
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which has been evaluated in the previous section, 



Gp{h,t2) = Gl^{h-h) + Gp{h,t2) . (5.12) 

Lepton propagators and self-energies have large equilibrium contributions dominated by 
gauge interaction, with small corrections (9(A^), 

^Lk(tl,t2) = 521(tl -t2) +55Lk(tl,t2) , (5.13) 

nk(ti, t2) = n'k''(^i - + 5nk(ti, t2) , (5.14) 

which include CP-violating source terms and washout terms. Clearly, inserting and 
5*^*^ in Eq. (|5.8I) must yield L^^(t,t) = 0, since no asymmetry is generated in thermal 
equilibrium!^ As discussed in Section 2, we also neglect washout terms for simplicity. One 
then obtains for the lepton number matrix Lk(t, t) to leading order in A, 

Lk(t,t) = z / [ dhtT[ST]^{h,t2)SX{t2-h) 
Jo Jo 

-6U^{t„t2)SZ-ih-h)] . (5.15) 



Here 5TTk is given by the two-loop graphs shown in Fig. [TJ which have to be evaluated 
with equilibrium propagators for lepton and Higgs fields and the nonequilibrium Majorana 
neutrino propagator. 

The equilibrium propagators with standard model gauge interactions remain to be 
evaluated. In the quasi-particle approximation one simply replaces energies k by complex 
quasi-particle energies = (k^ + m'^hY^'^ + i'y{k). In the following we shall consider 
two approximations: free equilibrium propagators with zero chemical potential as given in 
Eqs. dilD, (jA3) and (jMD, (DD, 

K^^iy) = ^iiy) , s^^iy) = s^M , (s-ie) 

and, as a rough approximation to full thermal propagators, free equilibrium propagators 
modified by thermal damping rates, 

K^iy) = A^(y)e-^*l^l , S^^iy) = S^,(y)e-^'l^l . (5.17) 

Remarkably, thermal widths turn out to be qualitatively more important than thermal 
masses, as we shall explain in Section 6. 

The two contributions to the self-energy 6V\kij (cf. Fig. [7]), 

5T^^^Atl,h) = ng(ti,t2) + ng(ti,t2) , (5.18) 



^°Note that thermal equihbrium does not correspond to a Gaussian state [42]. Therefore one has to 
include contributions from n-point functions which are not determined by equilibrium 2-point functions. 
However, such terms do not contribute to leading order in the Yukawa coupling A. 
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factorize into a product of Yukawa couplings, which contains the flavour dependence, and 
a trace of thermal propagators, 

nS(ti,t2) = -3^*1 (r^A*)^.,n«(ti,^2) , (5.19) 
ng(ti,t2) = S^ iv*X),, X,iT]S\h,t2) . (5.20) 

In the case of free equilibrium propagators for lepton and Higgs fields, we obtain for the 
self-energies V\^'^^^ and T\^'^^^: 



POO P 

nL'^>(ti,t2) = dh J 



dh J 



(27r)3 (27r)3 

- Gf{ti,h)S<{t2 - h)A<{t2 - h)A<{t2 - ti)]PL , (5.21) 
' d\ d^q' 



(27r)3 (27r)3 



dh j 



X [Gl\t,, h)S>{t2 - h)A>{t2 - h)Al{t2 - h) 

- G<{h,t,)Sl?{t2 - h)A^^{t2 - h)A>{t2 - ti)]PL , (5.22) 
d^q d^q' 
(27r)3 (27r)3 

X [G<it2, h)Sl^its - ti)A5(t3 - ti)A<(t2 - h) 

- Gl\t2,ts)S<{t, - ti)A< (t3 - h)A<{t2 - h)]PL , (5.23) 



poo p 

ni,'^<(ti,t2) = dh j 



(27r)3 (27r)3 

X [Gf{t2, t,)S>{t, - ti)A> (ts - ti)A>(t2 - ti) 

- G>(t2, h)Sl}{h - h)A\]{t, - ti)A>(t2 - h)]PL . (5.24) 

Due to the chiral projections at the vertices, only the scalar parts of the nonequilibrium 
Majorana propagators contribute, which are the same for G^, G^, G^ and G^ (cf. Eqs. 
flXT9|) - (1X241) ). 

PiGp(t, t')GPL = Gp(t, t')PL , Gp(t, t') = — cos(a;p(t - t'))r^\ujp)e-'^-'^'+''y' . 

UJp 

(5.25) 

The number of terms which contribute to the asymmetry L\^{t,t) can be significantly 
reduced by means of the following symmetry properties of the massless propagators: 



S^{yr = CS<iy)G-' , Si\yr = CSl\y)C-' , (5.26) 

A>(y)* = A<{y) , A\\yr = Afi-y) , (5.27) 
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Skiy) = 755-k(-l/)75 
A<(y)=A>(-y), 



sl\y) = i5S'\i-yH 

A\\y) = Alli-y). 



(5.28) 
(5.29) 



Employing these transformation properties one can derive the following useful relations 
among different contributions to the integrand of Eq. (I5.1ip : 



tr 



nt'^>(ti,t2)5<(t2-ti) 



tr 



-tr 
-tr 



n 



(1,2)< 



{tl,t2)S>{t2-h] 



ni?^^(^2,ti)5k>(ti-t2) 



(5.30) 
(5.31) 



Using these relations one obtains from Eq. ( 15. lip the compact expression for the lepton 
asymmetry 



LUt,t) = 12 lm{X*,{r]X*)J 

X dti dt2 Re ( tr 
Jo Jo 



(5.32) 



Since lm{X*^ (^-^*)ji} = 167rejj/(3M) (cf. Eq. fl2.12p ). the leading dependence of the flavour- 
diagonal lepton asymmetry L-\^ii{t, t) on the Yukawa couplings is identical to the dependence 
of the difference /Li(t, k) of lepton and anti- lepton distribution functions appearing in the 
Boltzmann equations. 

To proceed further in the evaluation of LkM(t, t), the following relation can be used to 
simplify the integrand, 

S'^{y)A'^{y)-S<{y)A<{y) = 



Qj-y) 

2q 



1 



l^coth 



^ j cos(%) cos{qy) — tanh ) sm{ky) sm{qy) 



.M 



k 



— I tanh I 4^ I cos{ky) sm(qy) + coth [ ^ ) sm{ky) cos{qy) 



. (5.33) 



One then obtains for the real part of the sum of products of thermal lepton and Higgs 
propagators (?/« 



Re(tr[(5^?(y3i)Aj(y3i) - S<{ys,)A<{y,,))S<{y,,)]A<{y,,)) 



Q{y 



13j 



16qq' 



coth (^Y^ (cos((A; + g)?/2i) + cos((A; - q)y2i] 
+ tanh (^cos((/c + g)?/2i) - cos((A; - q)y2i] 

(^coth (^^^ (cos((A;' + g')2/3i) + cos((fc' - gOz/si) 



+ tanh 



(3k' 



cos((A;' + g')z/3i) - cos((/c' - q')y: 
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+ -r-rT- I coth ( — ) ( sm((A; + g)?/2i) + sm((A; - q)y2i] 



kk' \ \2 

+ tanh (^sin((A; + g)?/2i) - sin((A; - q)y2i] 

(coth (^^^ (sin((A;' + g')y3i) + sin((A:' - q')y,,] 

+ tanh I — j (^sin((/c' + g')2/3i) - sin((A;' - q)yzi) 



(5.34) 



Defining the linear combinations of lepton and Higgs distribution functions (cf. (12.61) ). 



U^{k, q) = \- m + U{q) , fi^ik, q) = m + 
and using the relations 



coth ( ^ j + tanh 



coth 



one finds 



(5q 



tanh 



(3k 



2Mk,q) 



(5.35) 

(5.36) 
(5.37) 



Li,ii{t,t) = -en 32n f dti f dt2 [ dts f —f^{ujp)e~ 

Jo Jo Jo Vq,q' 



riti+ts) 



cos[ujpysi, 



X 



k - k' 

kk' 



fi^{k, q) cos{{k + q)y2i) + fi4,{k, q) cos{{k - q)y2i 
X ( fi^{k', q') cos{{k' + q')y23) + Mk', q') cos((fc' - q')y23) 
fi4,{k, q) sin((A; + g)?/2i) + //^(A;, q) sin((/c - q)y2i\ 

X (fi^{k\ q') sin((fc' + g>23) + h{k\ q') sin((A;' - g')l/23) 



(5.38) 



where we have again used the notation 
(27r)32g ' 

The functions fi^ and are well known from Weldon's analysis of discontinuities in 
finite-temperature field theory |35j. The sum of statistical factors 



fi^{k, q) = {l- m){l + U{q)) + mU{q) 



(5.39) 
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corresponds to decays and inverse decays of the massive Majorana neutrinos whereas 



Mk, q) = U{q)il - m) + mil + /^(g)) (5.40) 

accounts for their disappearance or appearance where a single quant, lepton or Higgs, is 
absorbed from or emitted into the thermal bath. The function fi^ contains the vacuum 
contribution, i.e., //<^ — > 1 as /3 — oo, whereas //<^ — > 0. 

We now have to perform the three time integrations in Eq. fl5.38p . It is convenient 
to express the products of cosine's and sine's as sum of products of exponentials. Each 
term then becomes a sum of four exponentials, where the energies u, k ± q and k' ± g' 
appear in different linear combinations, and the four complex conjugate exponentials. As 
an example, consider the integral 

X(t) = f dti [ dt2 [ ' dt3e-^^i*i+^^2t2+^n3t3g-E(ii+i3) ^ ^5 4^^) 

Jo Jo Jo 

with = Up — k — q, = Up — q' — k', and = fii— = k' + q' — k — q. A straightforward 
calculation yields 

-r fe"'^* + cos((]2t) - e"^ (cos(ilit) + cosin^t))) + 0{t) 
I{t) + T{t) = ^ ^ , (5.42) 

where 

0{t) = — ^ sin(n2t) - e"^ (sin(fiit) - smin^t)) ) (5.43) 

is of higher order in F at fii 3 = 0. Hence, this term does not contribute to the lepton 
asymmetry at leading order in F, i.e., in the Yukawa couplings. 

The two contributions in Eq. (15.381) . without and with the prefactor h-k.' /{kk'), add up 
to a single term proportional to k ■ k' / {kk') where k ■ k' denotes the product of 4-vectors. 
This is a consequence of Lorentz invariance of the vacuum contribution. The full result is 
now easily obtained from Eqs. f l5.38p and f l5.42p by adding the contributions with reversed 
sign of q and/or g', accompanied by the corresponding substitution /;$ — )■ Omitting 
the subleading terms O (cf. fl5.43p ). one finally obtains 

4 

LM,{t,t) = J2L^^j{t,t) , (5.44) 

a=l 

where 

L^t, t) = -eu Svr / /^^(a;p) ^F J] L^,^,^,(t; a, /3) (5.45) 
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and 



^k,q,q'(^;">/3) 



'^•^ (K-«(^ + ?)P + ?)(K-/3(A:' + g'))' + ¥) 

X (^e-^* + cos[(a(A; + q) - I3{k' + q'))t] 

- ( cos[(wp - a{k + q))t] + cos[(wp - (3{k' + g'))^])) , (5-46) 

^ ficp{k,q)fi4k',q') 

(K - «(A: - + ^)((a;p - + q')y + ^) 

X (^e^^* + cos[(a(A; - q) - /3{k' + q'))t] 

( cos[(wp - - + cos[(u;p - f3{k' + q'))t]) ) , (5.47) 

fi^{k,q)fi<p{k',q') 

(K - c^ik + q)f + f )(K - /3(A:' - g'))^ + f ) 
X (^e-^* + cos[(a(A; + q) - /3{k' - q'))t] 

- ( cos[(a;p - a{k + q))t] + cos[(a;p - I3{k' - q'))t]) ) , (5.48) 
f4 ^ m _ Mk,q)fi^{k',q') 

-^k,q,q'l'-i "5 PJ — 77 77 r2\// WrTi i\\2 



_ rt 

e 2 



^k,q.q'(^;«,/3) 



((cp - «(A: - q)y + - f3{k' - q')Y + ^^) 

X {e-^' + cos[(a(A; - g) - /3(A;' - g'))^] 

^ ( cos[(wp - a(A; - g))t] + cos[(u;p - P{k' - g'))t])) . (5.49) 



e 



This expression contains off-shell and memory effects which are not contained in Boltzmann 
equations. A detailed comparison will be given in the following section. 

So far we have neglected the thermal damping widths of lepton and Higgs fields due 
to gauge interactions, which are known to be much larger than the width of the heavy 
Majorana neutrino, 7/ ~ 7$ ~ g^T ^ A^M ~ F, for M < T. To estimate their effect we 
replace the free equilibrium propagators by 

Al^^iy) = A±(y)e-^*l^l , S^^^iy) = S^{y)e-^'^y\ . (5.50) 

This has a drastic effect on the calculation described above. For the dominant term in 
Eq. fl5.45p . -f'k^q^q' with a = (3 = 1, where the energy dominators can be C(F^), one now 
finds (7 = 7^ + 7$), 



Lkii{t,t) = -en 167r / 

J a 



k-k' 



q,q' kk'UJp 

17 

X 



((a;p - A; - g)2 + ^^){{ujp - k' - q'f + i^) 
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xfi<t>{k,q)Mk',q')f'^{ujp) 

xi(l-e-r*) , (5.51) 
where 7 = 7(/c, q) and 7' = '~f'{k', q'). Note that now all memory effects have disappeared. 

6 Boltzmann vs Kadanoff-Baym 

Let us now consider in detail the relation between the two results obtained for the lepton 
asymmetry: Eq. (12. lip from the Boltzmann equations and Eqs. (I5.44p - fl5.49p and (I5.5ip 
from the Kadanoff-Baym equations. 

Clearly, the overall CP asymmetry is identical in both cases and also the momentum 
integrations are very similar. Compared to the Boltzmann result the Kadanoff-Baym result 
has an additional statistical lepton-Higgs factor and expected off-shell energy denominators. 
Furthermore, there are 16 different terms corresponding to the various combinations of 
decay and inverse decay, appearance and dissappearance. The most striking difference 
is the time dependence of the integrand: the Boltzmann result has a simple exponential 
behaviour whereas the Kadanoff-Baym result has terms rapidly oscillating with time with 
frequencies 0{M) ^ F, a manifestation of memory effects. 

The time- dependence is contained in the integral X{t) given in Eq. fl5.4ip . Defining 

n, = Qi+'-r , = + , (e.i) 

and using the identities ^3 = ^1 + (^2 ^^i) + (^3 — ^2) and ^^2 = ^1 —^3, one has (cf. f l5.4ip ). 

I{t) = / rftie^^*^ / dt2i / rft32 e'^''^'+'^'''^ , (6.2) 
Jo J-h J-t2 

where tij = ti — tj. After performing the time-integrations, one obtains the result 
1 



(6.3) 



|^]l|2 V /V / n2ni 

which satisfies 

X(0) = X'(0) = X"(0) = , X'"(0) ^ . (6.4) 
For large times, t ^ 1/F, there remains a term oscillating with time. 



J_ + {e'^'' - 1) 

1 1^]^ I ^ ^^2^1 



(6.5) 
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This is in contrast to the Bohzmann result whose time-dependence is given by 



1 - e-^* 

XbW = , (6.6) 

with 

Xb(0) = 0, X^(0)7^0, (6.7) 

and ^sit) ~ 1/r = const for large times t ^ l/F. 

Where is the Boltzmann result hidden in the Kadanoff-Baym result, and in which limit 
is it recovered? To answer this question it is instructive to consider a modified integral 
X(t), where thermal damping rates 7 ~ 7' ~ g'^T are included, which affect the dependence 
on the time differences 1^2 ~ ^1 1 and — ^2! (cf. Fig. [7]), 

lit) = / dhe-^^' / dt2i / dh2 e^^i*2i'Tl*2il e*f^3t32-7'|t32| _ g) 
Jq J-h J -t2 

Compared to Eq. (I5.4ip the main difference is that the damping term in the t2i-integration 
changes sign at ^21 = 0. This is in contrast to the damping due to the Majorana neutrino 
decay width V. 

Carrying out the time-integrations one now obtains the result 



+ 1' 



1 g(ini-7)t /g(-jn*+7)t _ -j^ 



(zfii -7)(-zl]t +7) 
(zl]i+7)(-znt-7) V ) 



1 g(jn2-7-7')t _ 1 



(,fi2-7-y)(-^^i + 7) 

(zfi2 + 7-7')(-^^t -7) ^ ^ 

27 1 - e~^* 27 / (_jn*_y)i _ \ 

"n? + 72 r {ini + 7')((2fi2 - I'Y - 7') V ) 



(6.9) 

The first four terms reduce to Eq. (16.31) for 7 = 7' = 0. Particularly interesting is the 
last line in Eq. (16.91) . which is a contribution from the point ^21 = ^2 ~ = 0, where 
the damping term changes sign. This local contribution contains the only term which is 
enhanced by l/F and has Boltzmann-like time-dependence, 

lit) D Isit) = ^ ~ ^ . (6.10) 
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Note that as consequence of thermal damping all oscillatory terms are exponentially sup- 
pressed for times t > I/7, 



1 r 27 l-e"^* 27 




(6.11) 



{in, + 7) (-^nj - 7) (z^a + 7 - - 



7) 



The Boltzmann-like term fl6.10p . which originates from the point t2 = ti, vanishes for 



What is the order of magnitude of the lepton asymmetry (I5.44p relative to the Boltz- 
mann result in the case 7 = 7' = 0? The Kadanoff-Baym result depends on r = Ft, 
like the Boltzmann result, and in addition on the dimensionless parameter F/M ^ 1. In 
appendix C we shown that 



Hence, in this zero-width limit, due to rapid oscillations of the integrand, the lepton asym- 
metry obtained from the Kadanoff-Baym equation is at least 0{r/M) relative to the Boltz- 
mann lepton asymmetry. 

We are thus led to the conclusion that the lepton asymmetry obtained from the 
Kadanoff-Baym equations does not contain the Boltzmann result as limiting case as long 
as free equilibrium propagators are used for lepton and Higgs fields. This may not be too 
surprising. After all, the underlying assumption in our calculation has been that (gauge) 
interactions, much faster than heavy neutrino decay, establish kinetic equilibrium for lep- 
tons and Higgs particles. These interactions will unavoidably lead to thermal damping 
widths much larger than F. If these interactions are not taken into account in the calcula- 
tion of the lepton asymmetry, one misses the main contribution and obtains a misleading 
result. This means that at present the best estimate for the full quantum mechanical lep- 
ton asymmetry is given by Eq. ( I5.5ip . which leads to a temperature dependent suppression 
compared to the Boltzmann result. 

Note that the proposed incorporation of thermal damping rates leads to a Boltzmann- 
like result, Eq. (I5.5ip . which is valid for t > 1/F. For t < 1/F, all terms have to be kept, 
and one has dtLi^{t, t)\t=o = 0, which is a property of the exact result fIS.Sp . contrary to the 
Boltzmann approximation. 

7 Numerical analysis 

Let us now quantitatively compare the Boltzmann result 02.1 4p for the lepton asymmetry 



7 



0. 



^0 , 



for > , T = Tt fixed . 

M 



(6.12) 



fUt,k) = Ut,k)-fUt,k) 
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with the KadanofF-Baym result for the lepton asymmetry 

L^uit,t) = -tr[^oSt^,,it,t)] . (7.1) 

For free fields in thermal equilibrium both expressions are identical. For the Kadanoff- 

Baym result we use Eq. (I5.5ip which includes the estimated effect of thermal widths for 
lepton and Higgs fields. 

As shown in Appendix C, the Boltzmann result fl2.14p can be reduced to a two- 
dimensional momentum integral (cf. (IC.ISP ). 

fL^it, k) = FB(fc, /3) ^ (1 - e-^*) , (7.2) 
where we have defined 



1 /-OO /•'«maxW 1 

Fb(A;,/3) = - / dp k'dk'— 

X ^1 ) Mk, Up - k)f^ (cp) ; (7.3) 



here Up = a/ + , the bracket represents the product of 4-vectors divided by the 
corresponding energies, k ■ k'/{kk'), and the integration boundaries are 

n.\ _ -4:k'^\ u _ ^p-P y / ^ _ + P /y.N 

The dependence on temperature (/3 = 1/T) enters through the equilibrium distribution 
functions of Higgs particles and leptons, 

f^{k, g) = 1 - fi{k) + f^{q) , q = ujp-k , (7.5) 

^'^^) = ;^r^' ^^(^) = ;^7rT' = ^ ■ (7.6) 

The Kadanoff-Baym result f l5.5ip for the lepton asymmetry, which includes effects of 
thermal damping, takes the same form as the Boltzmann result 

LUt,t) = -^FMk,(3)l{l-e'^') . (7.7) 

Since the integrand of the momentum integrations contains two delta-functions less than 
the expression for the Boltzmann result, the function Fkb(^,/3) can only be written as a 
four-dimensional integral (cf. (1C.19P ). 



FKB{k,[5) = —- dp k'dk dq dq'— 
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fi^{k,q)fi^{k',q')f^{ujp) 



(7.8) 



(K - A; - g)2 + 72)((^p _ k' - q'f + y^) 



with the integration boundaries 



q± = \p±k\ , q'^ = \p±k' 



(7.9) 



For the thermal widths we use the estimate 7 ~ 7' ~ ~ 0.1 T (cf. |11]). Note that 
the damping in a non-Abehan plasma is considerably stronger than in an electromagnetic 
plasma at the same temperature. 

It is instructive to compare the Boltzmann and Kadanoff-Baym results with the pre- 
diction of quantum Boltzmann equations. As shown in [|15 | I2U ] . these equations lead to 
an additional statistical factor compared to Boltzmann equations, which implies for the 
lepton asymmetry 



In [T9||20|. this enhancement has been included in an effective, temperature-dependent CP 
asymmetry. 

In Fig. E] Boltzmann and Kadanoff-Baym results for the lepton asymmetry are com- 
pared. At momenta k ~ 0.2, where both distributions peak, the differences are less than 
20%, at larger momenta they reach at most 50% (cf. Fig. ITOl) . At temperatures T ~ 0.3, 
where leptogenesis takes place for typical neutrino parameters [30 ^ i45 j . differences are es- 
sentially negligible. 

Boltzmann and quantum Boltzmann results for the lepton asymmetry are compared 
in Fig. [9l At momenta k ~ 0.2, where both distributions are maximal, the differences 
can exceed 100%, and they remain large also at larger momenta (cf. FigJTOl). An enhance- 
ment (9(100%) at T ~ 1 is qualitatively consistent with the enhancement found for the 
temperature-dependent CP asymmetries in [T9l[20] . 

The Kadanoff-Baym result strongly depends on the size of the thermal damping rates. 
For 7, 7' — )■ 0, off-shell effects dissappear, and the Kadanoff-Baym result approaches the 
quantum Boltzmann result. Numerically, already for 7 ~ 7' ~ 0.01 T the differences 
are negligible. However, in a non-Abelian plasma, damping rates are large and, as a 
consequence, they almost compensate the enhancement due to the additional statistical 
factor contained in the quantum Boltzmann as well as the Kadanoff-Baym result. We 
conclude that, according to our estimates, the conventional Boltzmann equations provide 
rather accurate predictions for the lepton asymmetry. 




(7.10) 
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2ujpk - 2ujpk' - M2 
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fi<j,{k,Up - k)fi^{k',Up - k')f^{ujp) 
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Figure 8: Comparison of the lepton asymmetry distribution functions obtained from Boltz- 
mann equations (B, dot-dashed hne) and Kadanoff-Baym equations (KB, full line) for three 
different temperatures; temperature and momentum are given in units of M. 
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Figure 9: Comparison of the lepton asymmetry distribution functions obtained from Boltz- 
mann equations (B, dot-dashed hne) and quantum Boltzmann equations (QB, dashed hne) 
for three different temperatures; temperature and momentum are given in units of M. 
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Figure 10: Ratio of Kadanoff-Baym and Boltzmann lepton asymmetries (upper panel) and 
ratio of quantum Boltzmann and Boltzmann lepton asymmetries (lower panel) for three 
different temperatures; temperature and momentum are given in units of M. 
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8 Summary and conclusions 



The goal of leptogenesis is the prediction of the cosmological baryon asymmetry, given 
neutrino masses an mixings. In a 'theory of leptogenesis', it must be possible to quantify 
the theoretical error on this prediction. This requires to go beyond Boltzmann as well as 
quantum Boltzmann equations, such that the size of memory and off-shell effects can be 
systematically computed. 

In the present paper we have shown how to calculate the lepton asymmetry from first 
principles, i.e., in the framework of nonequilibrium quantum field theory. Our calculation 
is entirely based on Green's functions, and it therefore avoids all assumptions which are 
needed to arrive at Boltzmann equations. 

Two key ingredients make the problem solvable. First, the thermal bath has a large 
number of degrees of freedom, all standard model particles, compared to only one parti- 
cle out of equilibrium, the heavy neutrino. Hence, the backreaction of its equilibration 
on the temperature of the thermal bath can be neglected. Second, the heavy neutrino 
is only weakly coupled to the thermal bath and we can use perturbation theory in the 
corresponding Yukawa coupling A. 

The weak coupling of the heavy neutrino to the bath allowed us to obtain analytic 
expressions for the spectral function, which do not depend on initial conditions, and for the 
statistical propagator. In Section 4 we have discussed two solutions of the Kadanoff-Baym 
equations, which correspond to thermal and vacuum initial conditions. The statistical 
propagator which interpolates between vacuum at t = and thermal equilibrium at large 
times can then be used in the computation of the lepton asymmetry. 

Thermal leptogenesis has two vastly different scales, the width F of the heavy neutrino 
on one side, and its mass M, temperature T of the bath and thermal damping widths 7 
on the other side, 

F ~ A^M < 7 ~ ^^T < T < M . 

Typical leptogenesis parameters (cf. [2j) are F ~ 10^'' M, 7 ~ 0.1 T, T ~ 0.3 M, M ~ 
lO^'' GeV. The existence of interactions in the plasma, which are fast compared to the 
equilibration time rjy = 1/F of the heavy neutrino, is always implicitly assumed to justify 
the use of Boltzmann equations for the calculation of the asymmetry, but their effects are 
usually not explicitly taken into account. 

The main result of this paper is the computation of the lepton asymmetry in Section 5, 
where the nonequilibrium propagators of the heavy neutrino and free equilibrium propaga- 
tors for massless lepton and Higgs fields are used. Compared to Boltzmann and quantum 
Boltzmann equations, the crucial difference of the result (I5.44p - fl5.49p are the memory 
effects, oscillations with frequencies (9(M), much faster than the heavy neutrino equilibra- 
tion time tn = 1/F. These oscillations strongly suppress the generated lepton asymmetry 
Lk{t,t) compared to the Boltzmann result fi{t,k). In fact, as shown in appendix C, the 
ratio Lif_{t,t)/ fi^lt, k) vanishes in the 'zero-width' limit F/M — 0, with r = Ft fixed. 
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This situation changes when the interactions, which in the Boltzmann approach are 
assumed to estabhsh kinetic equihbrium, are exphcitly included in the calculation. Lepton 
and Higgs fields in the thermal bath then acquire large thermal damping widths 7 ~ g'^T, 
which cut off the oscillations. As a consequence, the predicted lepton asymmetry is similar 
to the quantum Boltzmann result, except for off-shell effects which are now included. For 
small damping widths, 7 <^ T, the off-shell effects are negligible. They are large, however, 
in the standard model plasma. According to our calculation, using 7 ~ 0.1 T, the damping 
effects essentially compensate the enhancement due to the additional statistical factor of 
the quantum Boltzmann equations. We conclude that, after all corrections are taken into 
account, the conventional Boltzmann equations again provide rather accurate predictions 
for the lepton asymmetry. Note that the classical Boltzmann behaviour emerges at large 
times, t > 1/r > 1/7, while at early times all terms are of similar magnitude, and all 
quantum effects have to be kept. 

As already emphasized in [8], it is of crucial importance to include gauge interactions in 
the Kadanoff-Baym approach to make further progress towards a 'theory of leptogenesis'. 
It remains to be seen whether the qualitative effects of thermal damping, as discussed in 
this paper, will then be confirmed or whether new surprises are encountered. 
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A Thermal propagators 

In the following we list all propagators, which are needed in the calculation described in 
Section 5, as functions of relative time y = ti — t2 and total time t = {ti + t2)/2. 



Free massive scalar (a;q = y/w? -\- q^) 

Aq(y) = — sin^y) , (A.l) 



K^y) = ^ coth ) cosKy) , (A.2) 



^1 



Ai^(2/) = ^ (^coth (^^^ cos^y) - .sinKli/l)^ (A.3) 



= A+(y)-%ign(y)A-(y) , 



(y) = ^ (^coth (^^^ cosKy) + zsinKlyl)^ (A.4) 



A+(t/) + ^sign(|/)Aq(7/) , 



^^^^-^ " 2^ (^°*^^ (^) ^^^('^'i^) - i sin(a;q7/)^ , (A.5) 
Aq(y) = ^ (^coth (^^^ cos(a;qy) + isin(a;qy)^ . (A.6) 

Free massive Dirac fermion {uj\^ — vp? + k^) 

ley 

'S'k (y) = i7ocos((x;ky) H sin(a;ky) , (A.7) 

'S'k (y) = -^tanh (^^^^ ^i7o sin(<^kZ/) - ^ J^'^ cos(c^kZ/)^ , (A.8) 
'S'k^(?/) = y ^cos(a;k|/)sign(|/) - itanh {^^i^ sin(a;k|/) 

+ "^^^^^ l^tanh {^^Y^ cos{uky) - i sin(a;k|y|) ) (A.9) 
= - ^sign(y)5k-(y) , 

Skiv) = y ^- cos ((^k|/) sign (y) - itanh (^^^^ sin(a;ky) 

H ^^^^ — I tanh {—^ 1 cos(a;k|/) + ism{ujk\y\) ] (A. 10) 
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+ 



To 
2 

m — kj 



S^{y) = y (^cos(cjky) - itanh 



sin(cjk|/) 



2a;k ( 



tanh 



cosiuku) - i sin(a;ky) ) , 



7o 



S^(y) = — — cos(a;ky) — itanh 



+ 



2 

m — k7 
2u;u 



/3wk\ . , . 
sm(a;ky) 



tanh 



cos{uJky) + isin(a;ky) ) . 



The propagators for a massless left-handed fermion are obtained by the 
cjk ^ = |k|, S^- ^ PlS^-, where Pl = (1 - 75)/2. 



Free massive Majorana fermion (a;p = ^/M^Tp^) 



Gp (y) = ho cos{ujpy) + 



M - 



sm hx", 



.y))c-' , 



G+(l/) = -^tanh(^ 



Gl\y) 



y ( cos(a;py)sign(y) - itanh 



ijosm{upy) - ^ ^ ^ cos{LJpy)^ C ^ 
sm{ujpy) 



I3up 



M-py f 



tanh 



I3ujp 



cos{ojpy) - i sin{u}p\y\) 



C 



-1 



Gf{y) 



To 
2 

M - p7 



— cos (a;py) sign (y) — i tanh 



P * sin(a;py) 



r 



tanh 



2a;p V 

icos{idpy) — itanh 



cos(a;py) + i sm{ujp\y\) 



C 



-1 



sm{ujpy) 



M - p7 



2a;p V 

7o 



( tanh 



cos{upy) -ism{upy) 



H ^ 

2cc;p V 



- cos(a;py) - itanh ^^-^^ sin(a;py) 



tanh 



cos(a;py) + i sin(a;pt/) 
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Nonequilibrium massive Majorana fermion (interpolation between vacuum at t = 
y = and thermal equilibrium at t = oo, and memory integral) 

G^{y) = (^lo cosKy) + ^^^^ sinKy)^ e-rpl.l/2c<-i , (A.19) 

Gp{t,y) = - ^27osin(wp|/) - ^^—^ cos{ujpy)^ 

X tanh e-^pl^l/^ + r^\cOp)e''^^^ C'' , (A.20) 

G'^{t,y) = G;{t,y) - '-sign{y)G-{y) , (A.21) 

Gf{t, y) = G;{t, y) + ^sign(y)Gp (y) , (A.22) 

G>(t,y) = G+(t,y)-^Gp(i/), (A.23) 

G<(t,y)=G+(t,y) + ^G+(y) , (A.24) 



Cp.mcml^'^/) = -^tanh [^Yj f^7osin(wpy) - '^'^ cos^uj^y) 



■'p 



X 



(e-rpli/l/2 _ e-rp*) c^' . (A.25) 



B Feynman rules 



For completeness, we list in the following the Feynman rules for the Standard Model 
Lagrangian with right-handed neutrinos given in Eq. (11. 2p : a,/3 are spinor indices and 
a,b, . . . are SU(2) indices. 



Majorana neutrino 
N 



Lepton doublet 
/ 

• «• ■ 

,a,a,i 



Gap{Xl,X2) 
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• Higgs doublet 



^ I ^ 

• Vertices 




iKl^ab{PR)al3 



iXii{CPL)pa^ab 



iViji^ac^bd + ^ad^bc){PRC)aP 
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C Zero-width limit 



In this section we consider the Kadanoff-Baym result for the lepton asymmetry normahsed 
to the Boltzmann result, L]^{t,t)/ fiit, k), in the zero-width limit as defined in Eq. (I6.12p . 
i.e., 

r 

— ^ , T = Tt fixed . 
M 

To this end we have to evaluate the corresponding momentum integral fl5.45p in this limit. 



C.l Boltzmann equation 

Consider first the Boltzmann result for the lepton asymmetry given in Eq. fl2.14p . 



fUt, k) = -e,,^ / k-k' {27c'f 6\k + q - p)(27r) + q' - p) 

v/q,p,q',k' 

x/.^(A:,g)/^>p)^(l-e-^*) . (C.l^ 



r 

The integration over q and q' can be performed using the (5-functions, which leads to 



xMA:,g)/^>p)^(l-e-^*) , (C.2) 

where q = |q| and q' = |q'|. The product of 4-vectors, k ■ k' = kk'{l — k ■ k'), depends 
on the angles between the different momenta. It is convenient to define the angles with 
respect to the momentum p: 6 = Z(k, p), 6' = Z(k', p) and (f' = Z(k^, k'^); here k_|_ and 
k'j_ are perpendicular to the vector p, i.e., k = ky + k^ and k' = ky + k'j_. In terms of these 

angles the unit vectors k and k' are given by (see Fig. IVVp 

(C.3) 





f cose \ 




f cos 6' \ 




sin 9 




sin 9' cos (p' 




K J 




\ sin 9' sin Lp' 1 



with k • k' = cos 9 cos 9' + sin 9 sin 9' cos 99'. We then obtain 

n POO p1 p27r -i 

hi{t,k) = --^ / d'p / k'''dk' / dcos9' / V 



X (1 — cos 9 cos 9' — sin 9 sin 9' cos ip') (C.4) 
X5{k + q- cup)6ik' + q'- Up)Mk, q)r^'M^ (l - 6"^*) . 

47 




Momentum conservation relates the energies q and q' to p, A;, k' and the angles 6 and 6', 

g= |p-k| = {p^ + k'^ ~2pkcosey/^ , (C.5) 
g' = |p - k'l = + k'^ - 2pkcos0'y/^ . (C.6) 

We can now make use of rotational invariance of the distribution function, 

fUt,k) = ^ [ dn^fut,k) . (C.7) 



An „ 

Changing variables, 

dq = ——dcosO , dq' = —- — dcosO' , (Ci 
q q' 

one arrives at 

r,u,>r'. n'-.>(. p^+k^-q^p^+k'^-q^\ 



h^it,k) = -^l J dp k'dk' j'' dq j'^ dq' [l- 



2pk 2pk' 



1+ 

x—5{k + q-u^)5{k' + q'-uj^)fi^{k,q)r^%uj^)Ul-e-^') , (C.9) 
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where the hmits of integration are given by the maximal and minimal value of q and g', 
respectively, 

g± = |fc±p| , q'^ = \k' ±p\ . (C.IO) 

Consider now the argument of one 5-function, VLi = Up — k — q, with fi™™ = ^p — k — q+ 
and ^5"^^ = Up-k-q^ (cf. Eg. (ICTD ) . Obviously, the conditions nf'' < and fi^^^^ > 
limit the integration range in p for given momentum k, 

P>- ^-=p^Uk). (C.ll) 

Similarly, the constraint p > (Af^ — Ak"^) / {Ak') restricts the integration range in k' for 
given p, 

k' > ^ k'^ Jp) , k' < ^ k'^M ■ (C.12) 



Changing again variables from q and q' to VLi and fia, respectively, and using 

d{p,k',n^,a,) _^ 

dip, k', g, q') 
the integral can now be written as 



(C.13) 



max 



fLiit,k) = -^- I dpi dk' I dn, I dQs 



mm 



Pmin(fc) •^k'miu'yP) 



ujp \ 2pk 2pk' 

xMfc,g)/^>p)^(l-e-r*) . (C.14) 

The limits of integration have been chosen such that they contain the points i7i = and 
fis = 0, which correspond to energy conservation, q = Up — k and q' = Up — k', respectively. 
Hence, the integration on Qi and Q3 can trivially be carried out, and we obtain the final 
result 

^ , e,,l r , 1 / 2upk-M^2ujpk' -M^ 
fUt, k) = - — - / dp dk —[1 — i^— 

X MKup - k)r^\oOp)^ (1 - e-^*) . (C.15) 
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C.2 Kadanoff-Baym equation 

We are now ready to evaluate the leading contribution of the Kadanoff-Baym result for 
the lepton asymmetry. It is given by Eq. (15.461) with a = P = 1, and it can be written in 
the form 



Lkiiit,t) = -en Stt 



k ■ k' 

ir 

2" — cos((u;p — k — j ^e~^ — cos(((X'p — k' — q')t) 



— sin((a;p — k — q)t) sin((a;p — k' — q')t) 



(C.16) 



We first change variables, (q, q') — )■ (p, k'), with p = q + k = q' + k', and use rotational 
invariance, 



1 r 

Lkiiitjt) = ^ / "^^k Lkii(t,t) . 
Choosing again angles according to Fig. \'r!\ the integral (1C.16P becomes 

/I /»1 /»27r 

dcos6 / dcosO' / d(p' (1 — cos 6* cos 6*' 
1 J-i Jo 

-sin^sin^'cos(^') F{d,9\---) 
= i2nfj dcosO j dcose' il-cosecose')F{e,e',---) . 



(C.17) 



(C.18) 



where we have used that the function F{6, 6',- ■ ■) does not depend on the angle As 
in the previous section, we now change the integration variables from {9, 6') to (g, q'), and 
using Eq. (IC.SP we obtain 



Liiii{t, k) 



en 1 



'•'=max(p) r<i+ ri+ 1 

dp I k'dk' / dq / dq'- 



X I 1 - P" + k\-^"p" + k'\-^'" \ f,,ik,q)Uk',q')r^^M 



2pk 



2pk' 



X 



2^- 
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X 



2* — cos(((X'p — k — q)t)J (^e " — cos((u;p — k' — q')t)J 



— sin((ci;p — k — q)t) sin((ci;p — k' — q')t) 



(C.19) 



where the hmits of integration are given in Eqs. (ICnoj) - (ICl2|) . We have restricted the 
integration over p and k to the range for which the intervals and [qL^q'^ contain 

points satisfying Wp — /c — g = and Wp — /c' — g' = 0, respectively. This finite part of the 
integral could then be (9(l/r), which is required to match the Boltzmann result for the 
lepton asymmetry. The remaining part is 0{1) and therefore suppressed compared to the 
Boltzmann result. 

Remarkably, the integral (lU.lOp is a sum of terms each of which factorizes into a product 
where one factor depends on q but not on g', whereas the other factor depends on q' but 
not on q. Hence one obtains 



L^ii{t,t)(x j dp k'dk' y2Vi{q^,q+)Qi{q'_,q'_^) 



(C.20) 



where we have dropped the dependence of the factors Vi and Qi on k, p and k' for simplicity. 
Because of the factorization, we can now perform the integrations on q and q' separately. 

Naively, one may think that in the zero-width limit F/M — )■ the cosine terms can 
be set to one. But for large time t, they oscillate fast, which leads to a different result. 
Consider the following contribution to the integral flaTOj) . 



V{q-,q^ 



1+ 



dq- 



F{q) 



^p-k-qy + — 



^2 cos((wp -k- q)t) 



where F{q) has no poles. Changing the integration variable from g to ^ 
= OJ-p — k — q, one obtains 



(C.21) 
2Vti/T, with 



'P (^^mm: 2^max) ^^j^ / ( '^P 



k Z 

2 



X 



z + 1 



(C.22) 



where Zmin = 2fi™™/r and 2;max = 2Q^^^/r, with z^i^ < and ^max > 0. In the limit 
r/M — > with T = rt fixed, the integration limits approach Zmin — >■ — c>o and 2;max — ^ +oo, 
respectively. The integral is now easily evaluated by means of the residue theorem leading 
to the result 



rviz, 



rami ^max^ 



-TT 



r-s>o 



F u- 



k — i— 



+ F {cj^-k + i- 



r-^o 
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= -27TF{ujp - k)e'^ 



(C.23) 



In Eq. flal9|) the term V appears together with a second term, 



r<i+ 

V'{q-,q+)= / 

J (7_ 



dq- 



F{q) 



UJp-k- qf + ^■ 

which can be evaluated in the same way as V in the zero- width hmit, yielding 



(C.24) 



TV'iz, 



miri) ^max ) 



IT 



r^o 



k - i- ] + F { u- 



, .r 

k + t- 
2 



e 2 



r^o 



2nF{ujp - k) e 



(C.25) 



Clearly, the two terms V and V add up to zero. The same result is obtained for the second 
factor Q after the q' integration, as well as for the product of two sinus functions. 

We conclude that the integral (1C.16I) does not contain a contribution 0{1/T). Hence, 
the ratio of Kadanoff-Baym result and Boltzmann result, Lk(t, t) / fiit, k), approaches zero 
in the limit T/M -> 0, r = Ft fixed. 



D Equilibrium contribution 

In Section 5 we argued that the equilibrium part of the heavy neutrino propagator does 
not contribute to the lepton asymmetry. In this section we verify this claim. 

The heavy neutrino propagator has an equilibrium and a nonequilibrium part, 

Gp{h,h) = G'^{h-t^) + Gp{h,t3) , (D.l) 

whose main difference lies in the time dependence, 

G'^ih - h) oc , G'p(ti, h) oc e-5(*^+*^) . (D.2) 

The computation of the lepton asymmetry in Section 5 was based on the nonequilbrium 
part, and it involved the time integral X (cf. Eq. (I5.4ip ). Because of the different time 
dependence given in Eq. (ID. 20 . the contribution of the equilibrium part to the asymmetry 
involves instead the integral 

J{t) = [ dti [ dt2 [ ' dh e-*^i*i+*^2t2+i03i3g-||<i-t3l ^ p 3) 

Jo Jo Jo 

which differs from X only with respect to the damping factor. Qi, Q2 and Q3 are different 
linear combinations of energies, which satisfy Qi = Q2 + ^3- 

In order to evaluate the integral JT", we have to split the time integration. 



J{t) = [ dti r dt2 r dh e-^^i*^+*^2*2+*^3*3e-^(*^-*^) 
Jo L Jo Jo 
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+ 



+ 



t2 



dU e 



-jS7ltl+jr22<2+if^3i3 -E(t3-ti) 



[D.4) 



Note the change of sign in the damping factor of the last two terms. As in Section 5, it 
is convenient to use the variables = i7i — |r and ^3 = — |r, for which the integral 
simplifies to 



Jit) 



dti 



+ 



dU 



dU 



t2 



tl 



g-if2iti+if22t2+«f^3i3 



^^^g-*f2ltl+«f22t2+«S^3*3 _j_ 



[D.5) 



t2 



Performing the ^3 integral and using the relation Qi = Q2 + ^3, we obtain 



Jit) 



dti 



+ / dt2{e 



1 

1 



Jft2t2 



_l_ -in*t-i ^ ( iQ.lt2 _ iO^ii 2^2*2 ^ 

^ e e J 



1 e 



i^2t2 



[D.6) 



It is now straightforward to carry out the integrations over ti and which leads to 

2 



J it) + J* it) 



V (fii + ^3) {cos^iVti - VL^)t) - 1 + {cosiVLit) - cos(fi3t))e~^ 



D.7) 



r2 



sin((^7l — ^^3)^) — (sin(i7it) — sin(i73t))e^ 



Note that the expression has no pole at i7i = ^3. 

As in appendix C we now have to evaluate the momentum integral 



S 



dVLi 



dns iJ + T) 



(D.8) 



with the integration limits given below Eq. (IC.lOp . To perform the zero-width limit, we 
again introduce the variables ^1,3 = 2fii^3/r. For F — )■ 0, the limits of integration 2™™ and 
approach — cxo and +00, respectively. The 2;3-integration can now be carried out by 
means of the residue theorem. The integrand of the remaining Zi-integration has a double 
pole. The integration can again be performed using the residue theorem, and we find that 
TS approaches zero in the limit F/M — 0, r = Ft fixed. Hence, the equilibrium part of 
the heavy neutrino propagator does not contribute at leading order in F/M. 
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